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THE COMPLEX ZEROS OF RANDOM SUMS 
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This paper is dedicated to the memory of Larry Shepp. 


This paper extends earlier work on the distribution in the com¬ 
plex plane of the roots of random polynomials. In this paper, the ran¬ 
dom polynomials are generalized to random finite sums of given “ba¬ 
sis” functions. The basis functions are assumed to be entire functions 
that are real-valued on the real line. The coefficients are assumed 
to be independent identically distributed Normal (0,1) random vari¬ 
ables. An explicit formula for the density function is given in terms of 
the set of basis functions. We also consider some practical examples 
including Fourier series. In some cases, we derive an explicit formula 
for the limiting density as the number of terms in the sum tends to 
infinity. 


1. Introduction.. The problem of characterizing the distribution of 
the roots of random polynomials has a long history. In 1943, Kac [15] stud¬ 
ied the real roots of random polynomials with iid normal coefficients. He 
obtained an explicit formula for the density function for the distribution of 
the real roots. 

Following the initial work of Kac, a large body of research on zeros of 
random polynomials has appeared - see [2] for a fairly complete account of 
the early work in this area including an extensive list of references. Most 
of this early work focused on the real zeros; [5], [8] and [25] being a few 
notable exceptions. The paper of Edelman and Kostlan [4] gives a very 
elegent geometric treatment of the problem. 

In more recent years, the work has branched off in a number of directions. 
For example, in 1995, Larry Shepp and I derived an explicit formula for the 
distribution of the roots in the complex plane (see [22]) when the coefficients 
are assumed to be iid normal random variables. A short time later, Ibragimov 
and Zeitouni [11] took a different approach and were able to rederive our 
results and also find limiting distributions as the degree n tends to infinity 
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under more general distributional assumptions. See also [13] and [14]. 

Also in the late 1990’s, it was pointed out that understanding deeper sta¬ 
tistical properties of the random roots, such as Appoint correlations among 
the roots, was both interesting mathematically and had important implica¬ 
tions in physics (see, e.g., [20], [7] and [21]). 

A number of papers have appeared that attempt to prove certain spe¬ 
cific properties under increasingly general distributional assumptions. For 
example, in 2002, Dernbo et al. [3] derived a formula for the probability 
that none of the roots are real (assuming n is even, of course) in the case 
when the coefficients of the polynomial are iid but not necessarily normal. 
Other papers have continued to study real roots—see, e.g., [12]. Another 
property that has been actively studied in recent years is the fact that as 
n gets large the complex roots tend to distribute themselves close to and 
uniformly about the unit circle in the complex planes-see, e.g., the papers 
by Shiffman and Zelditch [23], Hughes and Nikeghbali [9], Ibragimov and 
Zaporozhets [10], Pritsker and Yeager [19] and Pritsker [18]. Also, Li and 
Wei [17] have considered harmonic polynomials-polynomials in the complex 
variable z and it’s conjugate z. 

Using a very different approach, Feldheim [6] has derived a result that 
with some work can be shown to be equivalent to the results presented 
herein. 

Recently, Tao and Vu [24], drawing on the close connection with random 
matrix theory, derived asymptotic formulas for the correlation functions of 
the roots of random polynomials. They specifically address the question of 
how many zeros are real. 

The results summarized above mostly establish certain properties of the 
roots under very general distributional assumptions. The price paid for that 
generality is that most results only hold asymptotically as n —> oo. In con¬ 
trast, this paper introduces a modest generalization to the core assumptions 
underlying the results in [22] and we show that analogous explicit formulas 
can still be derived for any value of n. Specifically, instead of considering 
polynomials, ]Cj=o r lj z3 ■ we generalize the “basis” functions z J to be any set 
of entire functions, fj(z ), that are real on the real line. So, to that end, we 
let 

n 

p n(z) = J2 r ljfji z ), zeC, 

3=0 

where n is a fixed integer, the r/j' s are independent identically distributed 
IV(0,1) random variables, and the functions fj are given entire functions 
that are real-valued on the real line. We derive an explicit formula for the 
expected number of zeros in any measurable subset H of the complex plane 
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C. The formula will be expressed in terms of the following functions 

A)(z) 

n 

= B o( z ) 

n 

= £|/j(*0I 2 > 


3 =0 

3=0 

Ai(z) 

TL 

= B i( z ) 

Tl 

= Y,f3( Z )fj( Z )’ 


3=0 

3=0 

A 2 (z ) 

Tl 

= ^2fj( z ) 2 , B 2 {z) 

TL 

= 


3=0 


and 



(1) 

Dq(z) = sj B 0 (z ) 2 - 


and, lastly, 




E\(z) = \J A-2 (z)Aq (z) — Ai(z) 2 . 


Notes: (1) As usual, an overbar denotes complex conjugation and primes 
denote differentiation with respect to z. (2) The function E\ will only be 
needed in places where the argument of the square root is a positive real. At 
such places, the square root is assumed to be a positive real. (3) Throughout 
the paper, we follow the usual convention of denoting the real and imaginary 
parts of a complex variable z by x and y. respectively, i.e., z = x + iy. 

Let v n (Si) denote the (random) number of zeros of P n in a set 12 in the 
complex plane. Our first theorem asserts that throughout most of the plane 
this random variable has a density with respect to Lebesgue measure on the 
plane: 


Theorem 1.1. For each measurable set 12 C {z £ C | Dq(z) / 0}, 

(2) Ei/ n (12) = / h n (x,y)dxdy , 

Jn 

where 

B 2 D'q ~ -B 0 (|L?i | 2 + |Ail 2 ) + (ApgiAi + 
kD'I 

It is easy to see that the density function h n is real valued. It is less 
obvious that it is nonnegative. We leave this sanity check to the reader. As 
we see from the above theorem, places where Dq vanishes are special and 
must be studied separately. The real axis is one such place: 


h n (z ) 
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Theorem 1.2. On the real line, the function Dq vanishes. For each 
measurable set 0 C IR, 



where 



In the case where the fj^zfi s are just powers of z, our results reduce to 


those given in [22] . The proof here parallels the analogous proof given in [22] 
but there are a few differences, the most important one being the fact that, 
in general, the function B\ is not real-valued like it was in [22]. It seems 
that this explicit formula has not been derived before and, as shown in later 
sections, there are interesting new examples that can now be solved. 

While the definition of h n in Theorem 1.1 looks rather complicated, it is 
nevertheless amenable both to computation and, with some choices of the 
functions fj , it is amenable to asymptotic analysis as well. 

In the following section, we derive the explicit formulas given above for the 
intensity functions h n and g n . Then, in Section 3, we look at some specific 
examples and finally, in Section 4, we offer some speculation and suggest 
future research directions. 

2. The Intensity Functions h n and g n .. This section is devoted to 
the proof of Theorems 1.1 and 1.2. We begin with the following proposition. 

Proposition 2.1. For each region E C whose boundary intersects the 
set {z | Dq(z) = 0} at most only finitely many times, 



( 4 ) 


where 


BiDq + B 0 B\ — A 0 Aj 

BqF>o + Bq — AqAq 


( 5 ) 


Proof. The argument principle (see, e.g., [1], p. 151) gives an explicit 
formula for the random variable namely 
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Taking expectations in (6) and then interchanging expectation and con¬ 
tour integration (the justification of which is tedious but doable), we get 


( 7 ) 


E^n(ft) = 1 . f 
2m Jd 


E^f\dz. 

an P n {z) 


The following Lemma shows that, away from the set {z \ D$(z) = 0}, the 
function 


( 8 ) 


F{z) = E 


Pn{z) 


simplifies to the expression given in (5) and, since we’ve assumed that dLl 
intersects this set at only finitely many points, this finishes the proof. J3 


Lemma 2.2. Let F denote the function defined by (8). For z 0 {z \ Dq(z) 


0}, 


F = 


BiD 0 + B 0 B i - A 0 Ai 


BqDo + Bq — AqAq 


Proof. Note that P n (z ) and P' n {z) are complex Gaussian random vari¬ 
ables. It is convenient to work with their real and imaginary parts, 

Pn(z) = £l + *&2» 

Pn( z ) = ?3 + 

which are just linear combinations of the original standard normal random 
variables: 

n n 

6 = Y, a Fh, 6 = V i>j i/j . 

j =0 j =0 

n n 

6 = 'E'CjVj, & = ^2 d jVj- 

3 =0 j =0 

The coefficients in these linear combinations are given by 
aj = R e(fj(z)) 
bj = Im (fj{z)) 

Cj = Re(/j(z)) 
dj = Im (/<(*)) 


fj ( z ) + /jO) 

2 

fj(z) ~ /jO) 

2 i _’ 

fj(z) + f'Az) 


fj(z) - f'Az) 


2 i 


( 9 ) 
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Put £ = [ £1 £2 £3 £4 ] T ■ The covariance among these four Gaussian 

random variables is easy to compute: 


( 10 ) 


Cov(£) = E££ t 


a T a 

a T b 

a T c 

a T d 

b T a 

b T b 

b T c 

b T d 

c T a 

c T b 

c T c 

c T d 

d T a 

d T b 

d T c 

d T d 


We now represent these four correlated Gaussian random variables in terms 
of four independent standard normals. To this end, we seek a lower trian¬ 
gular matrix L = [ Uj ] such that the vector £ is equal in distribution to 
L£, where C = [ Cl C2 C3 £4 Y is a vector of four independent standard 
normal random variables. The following simple calculation shows that L is 
the Cholesky factor for the covariance matrix: 

(11) Cov(£) = E££ t = E L(( t L t = LL T . 


Now, since £ = L£ and L is lower triangular (the symbol = denotes equality 
in distribution), we get that 

P{i( z ) _ £3 + ^£4 
Pn{z) £l+*£2 

D (hi + ^4l)Cl + (^32 + ^ 42 X 2 + (^33 + ^ 43 ) £3 + ^44 £4 

(In + 2/21X1 + H22C2 

Hence, exploiting the independence of the £j’s, we see that 


( 12 ) 


F(z) = E 


K(z) 

Pn(z) 


a£i + (3(2 
7£i + X2 ’ 


where 


a = Z31 + U41 fd = Z32 + 2/42 

7 = Zn + ihi S = ^22- 


Splitting up the numerator in (12) and exploiting the exchangeability of £1 
and £ 2 , we can rewrite the expectation as follows: 


F(z) = a fh/S) + ^ f(5 / 7 ), 

0 7 

where / is a complex-valued function defined on C \ IR by 


f(w) = E 


£1 

wCi + £2 ' 
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The expectation appearing in the definition of / can be explicitly computed. 
Indeed, 


fH = 


t 

l 


r 2 ir 


r 2 ir 


p cos 6 


I o wp cos 6 + p sin 6 
dQ 


o-P 


2 pdpdd 


w + tan 6 

and this last integral can be computed explicitly giving us 


/H 


1 

w + i ’ 

1 

w — i 


Im(ui) > 0, 
Im(itf) < 0. 


Recalling the definition of 5 and 7 , we see that 

7 _ hi _ Mi 
5 I22 I22 

In general, l\\ and I 22 are just nonnegative. However, it is not hard to show 
that they are both strictly positive whenever 2 has a nonzero imaginary part. 
Hence, 7 /5 lies in the lower half-plane, 5/7 lies in the upper half-plane, and 


(13) 


F(z) 


a 1 (5 1 

+ ^ l + i 

ia + fd 

*7 + 5 

h 2 — l a + + ^ 42 ) 

—hi + i{hi + I 22 ) 


At this point, we need explicit formulas for the elements of the Cholesky 
factor L. From (10) and (11), we see that 


a T a 

/ 2 

— Hi 



b T a 

= ^21^11 

b T b = 

7 2 1 l2 

l 21 + l 22 

c T a 

= ^ 31^11 

c T b = 

hihi + I 32 I 22 

d T a 

= Uihi 

d T b = 

Uihi + U 2 I 22 
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Solving these equations in succession, we get 


^11 = 

hi = 

a T a 

V a T a 

b T a 

h2 = 

(a T a){b T b) — (b T a) 2 

y/ a T a 

yja T aR 

hi = 

c T a 

h'i = 

(■ a T a)(c T b ) — (c r a)(6 T a) 

y/ a T a 

V a T aR 

hi = 

d T a 

I 42 = 

(■ a T a)(d T b ) — ( d T a){b T a) 

yj a T a 

y/ a T aR 

where 


R = \J ( a T a)(b T b ) 

— ( b T a ) 2 . 


Substituting these expressions into (13) and simplifying, we see that 
(14) 

—d T a + ic T a — i [a T a(—d T b + ic T b) — (— d T a + ic T a)b T aj /R 

^ ^ —b T a + ia T a + iR 

Recalling the definitions of Uj. bj, Cj, and dj given in (9), it is easy to check 
that the following identities hold: 


a T a = 

j(^o + 2i?o + ^o)> 



b T a = 

~l(A 0 - A 0 ), 

b T b = 

— \{Aq — 2Ro + ^4o), 

c T a = 

\{Ai + B\ + B\ + A\), 

c T b = 

-K^i-Rr + Ri-Ir) 

d T a = 

-{(A 1 + B 1 -B 1 -A 1 ), 

d T b = 

-\{Ai -B 1 -B 1 + A 1 ) 


Plugging these expressions into (14) and simplifying, we get that 

Ai + B\ + (A 0 Bi + B 0 Bi — AiB 0 — AqA\)/Dq 
Aq + B 0 + Dq 

where Dq is as given in (1). It turns out that further simplification occurs 
if we make the denominator real by the usual technique of multiplying and 
dividing by its complex conjugate. We leave out the algebraic details except 
to mention that a factor of Aq + 2Bq + Aq cancels out from the numerator 
and denominator leaving us with 

BiDq + BqB\ - A 0 Ai 
Dq(Bq + Dq) 



( 16 ) 
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or, expanding out Dq , 


(17) 


BiDq + B 0 Bi — A 0 Aj 
BqDq + Bq — AqAq 


□ 


Lemma 2.3. 
each a € IR, 


and 


On the real axis, F has a jump discontinuity. Indeed, for 

lim 

z—>-a: lm(.z)>0 Bq(^CL) 

lim F = 

z^a: Im(2:)<0 Bq^CL) 


Proof. Consider a point a on the real axis. On the reals, Ak = B \ c , for 
k = 0,1, and so Dq = 0. Hence, the right-hand side in (16) is an indetermi¬ 
nate form. To analyze the limiting behavior of F near the real axis, we first 
divide the numerator and denominator by Dq: 


B i + 


BqB i - A)Hi 


(18) 


F = 


D, 


'o 


Bq + Dq 


Now, only the ratio in the numerator is indeterminate. To study it, we start 
by expressing things in terms of the fj functions: 

BqBi - AoAi = J2 fj( z )fk ( z ) (, fj( z )fk(z ) - fj( z )fk(z ) 
j,k 

= fj( z )fk( z ) lm ( fj(z)fk(z )) 

j,k 


and 


D o = B$-\A 0 \ 2 = J2fj( z )fk( z )(fj( z )fk{ z )-fj(z)fk(z)) 

j,k 


fj( z )fk{z)lm ( fj( z )fk( z ] 

j,k 


Next, we write the first few terms of the Taylor series expansion of the / 7 ’s 
about the point z = a, substitute the expansions into the formulas above 
and then drop “high” order terms to derive the first few terms of the Taylor 
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expansions for BoB\ — AoA\ and Bq — |Ao| 2 . For the first expression, we 
only need to go to linear terms to get 

BqBx-AqAi = 2 iY^f j {a)fk{a){f'j{a)f k {a)-f j {a)f k {a)^y + o(z-a) 


= 2i (^Ai(a) 2 - A 0 (a)A 2 (a)) y + o(z — a) 

(as usual, we use y to denote the imaginary part of z). For the second 
expression, we need to go to quadratic terms. The result is 

D l = (/j( a ) 2 /fc(a) 2 - V 2 + °(( z ~ a ?) 

j,k 




^/j(a)/j(a)J I y 2 + o((z-a 


= 4 ^A 2 (a)A 0 (a) - ^i(«) 2 ) V 2 + °{{z — a) 2 ). 

Hence, we see that 

(19) B ° Bl ^° Al = -i E\(a) sgn(a) sgn(y) + o(z - a). 

Bo 

Combining (18) and (19), we get the desired limits expressing the jump 
discontinuity on the real axis. □ 


Proof of Theorem 1.1.. Without loss of generality, it suffices to consider 
regions H that are either regions that do not intersect the real axis or small 
rectangles centered on the real axis. We begin by considering a region H that 
does not intersect the real axis. Applying Stokes’ theorem to the expression 
for E/y n (fl) given in Proposition 2.1, we see that 

1 /• d 

Eu n (Q) = - / —F(z,z)dxdy. 
vr Jn dz 

Note that we are now writing F(z, z) to emphasize the fact that F depends 
on both z and z. Letting the dagger symbol stand for the derivative with 
respect to z, we see from Lemma 2.2 that 

( 20 ) 

— = {{BoDo + Bl-AoAo^BlDo + BxDl + fljSi + B 0 b\ - A^Ai) 

— {B\Do + BqBi — AoAi)(HqDo + -Bo Dq + 2BoBq — AgAo)| 
/(BoDo + Bl-\Ao\ 2 ) 2 . 
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Recall that we have assumed that the functions fj are entire and are real¬ 
valued on the real line. Hence, they have the property that fj{z) = fj(z). 
Their derivatives also have this property. Exploiting these facts, it is easy 
to check that 


( 21 ) 


Bl — B±, Aq — 2Ai, b\ — 1 ? 2 - 


Recalling that Dq = \JB^ — |^4q| 2 5 we get that 


( 22 ) 


d\ = 


BqB\ — A 0 Ai 

Do 


As explained in [22], substituting these formulas for the derivatives into the 
expression given above for dF/dz followed by careful algebraic simplifica¬ 
tions (see the appendix for the details) eventually leads to the fact that 
(I/ 7 i)dF(z,z)/dz equals the expression given for h n in the statement of the 
theorem. □ 

Proof of Theorem 1.2.. Consider a narrow rectangle that straddles an 
interval of the real axis: H = [a, b\ x [—e, e] where a < b and e > 0. Writing 
the contour integral for Er , n (H) given by Proposition 2.1 and letting e tend 
to 0, we see that 

Ez ' n ((a, b)) = —[ (F(x~) - F(x+)) dx, 

where i/ n ((a, b )) denotes the number of zeros in the interval (a, b) of the real 
axis and 

F(x-) = lim F(z) and F(x+) = lim F(z). 

From Lemma 2.3, we see that 


9 „W = A (F(l -)_ F(l+ )) = |h) 


This completes the proof. 


□ 


3 . Examples.. In this section, we consider some examples. The sim¬ 
plest example corresponds to the fj simply being the power functions: 

fj(z ) = z 3 . 






12 


ROBERT J. VANDERBEI 




Fig 1 . Random degree 10 polynomial: r/o + r/iz + r/ 2 Z 2 + ■ • ■ + gioz 10 . In this figure and the 
following ones, the left-hand plot is a grey-scale image of the intensity functions h„ and 
g n (the latter being concentrated on the x-axis). The right-hand plot shows 200,000 roots 
from randomly generated polynomials. Note that, for the left-hand plots, the grey-scales for 
h„ and g n are scaled separately and in such a way that both use the full range from white 
to black. 


As this case was studied carefully in [22], other that showing a particular 
example (n = 10) in Figure 1, we refer the reader to that previous paper for 
more information about this example. 

Each figure in this section shows two plots. On the left is a grey-scale 
plot of the intensity functions g n and h n . On the right is a plot of hundreds 
of thousands of zeros obtained by generating random sums and explicitly 
finding their zeros. 

The intensity plots appearing on the left were produced by partitioning 
the given square domain into a 440 by 440 grid of “pixels” and computing the 
intensity function in the center of each pixel. The grey-scale was computed 
by assigning black to the pixel with the smallest value and white to the pixel 
with the largest value and then linearly interpolating all values in between. 
This grey-scale computation was performed separately for h n and for g n 
(which appears only on the x-axis) and so no conclusions should be drawn 
comparing the intensity shown on the x-axis with that shown off from it. 
The applet used to produce these figures can be found at 

http://www.princeton.edu/~rvdb/JAVA/Roots/Roots.html 

Of course, the intensity function g n is one-dimensional and therefore it 
would be natural (and more informative) to make separate plots of values 
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of g n verses x , but such plots appear in many places (see, e.g., [16]) and so 
it seemed unnecessary to produce them here. 


3.1. Weyl Polynomials. Sums in which the ffl s are given by 


fj ( z ) = 


z 


3 


vn 


are called Weyl polynomials (also sometimes called flat polynomials). Figure 
4 shows the empirical distribution for the case where n = 10. For this case, 
the limiting forms of the various functions defining the densities are easy to 
compute: 


nm /io izj 

71 — 1-00 

lim A\(z) 

71—>00 

lim A 2 (z) 

71 — 1-00 

lim Dq(z) 

n—too 


= e 


= ze z 


lim B 0 (z) 

71—>00 


= e l-l 2 


ze\ z I 2 


(z 2 + l)e z2 

^f g 2(x 2 +y 2 ) _ g 4 (x 2 -y 2 ) 


lim B\(z) 

71—>00 

lim B- 2 {z) = (Izl 2 + l)e 

n—>oo 

lim E\(z) = e z 


Id 2 


The random Weyl polynomials are interesting because in the limit as n —> oo, 
the distribution of the real roots becomes uniform over the real line: 


Theorem 3.1. If fflz) = z- 7 /\ffl. for all j, then 


lim g n (x ) 

n—> oo 


1 

7r 


Proof. Follows trivially from Theorem 1.2 and the formulas above. D 


It is interesting to note that, in addition to the asymptotic uniformity 
of the distribution of the real roots, the complex roots are also much more 
uniformly distributed than was the case when we did not have the 1 j\ffl. 
factors. 

Figure 3 shows plots of g n for all of the examples considered here. 

3.2. Taylor Polynomials. Another obvious set of polynomials to consider 
are the random Taylor polynomials; i.e., those polynomials with 

«■> - f 

Figure 2 shows the n = 10 empirical distribution for these polynomials. 
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Fig 2. Random degree 10 Weyl polynomials: 170 + 1712 + 172 +• • ' + i? 80 ^ 7 g=j. T/ie empirical 

distribution on the right was generated using 2,500 random sums. 


3.3. Root-Binomial Polynomials. Let 



Figure 5 shows the n = 10 empirical distribution for these polynomials. This 
example is interesting because the real and complex density functions take 
on a rather simple explicit form. Indeed, it is easy to check that 


A 0 (z) 

= (1 + z 2 ) n , 

B 0 (z) 

= (i + N 2 ) n , 

Ai(z) 

= nz(l + z 2 ) n_1 , 

Bi(z) 

= nz( 1 + \z\ 2 ) n ~ 1 , 

M z ) 

= n(l + nz 2 )(l + z 2 ) n ~ 2 , 

B 2 (z) 

= n(l + n\z\ 2 ){\ + \z\ 2 ) n ~ 2 


The formula for the density on the real axis simplifies nicely: 


g{x) 


yfn 1 
7T 1 + X 2 


From this formula, we see that the expected number of real roots is y/n and 
that each real root has a Cauchy distribution. 
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Fig 3. The function g n for n = 10 for several choices of the fj’s 


3.4. Fourier Cosine Series. Now let’s consider a family of random sums 
that are not polynomials, namely, random (truncated) Fourier cosine series: 

fj(z ) = cos (jz). 

This case is interesting because these functions are real-valued not only on 
the real axis but on the imaginary axis as well: cos (iy) = cosh(y). Hence, Dq 
vanishes on both the real and the imaginary axes and, therefore, both axes 
have a density of zeros. The set of imaginary roots for a particular sum using 
the f/s map to a set of real roots if fj(z) is replaced with fj(z) = fj(iz ) = 
cosh( 2 :). Hence, the formula for the density on the imaginary axis is easy 
to compute by this simple rotation. The resulting density on the imaginary 
axis has this simple form: 

E/i 

nB 0 (iy )' 

An example with n = 10 is shown in Figure 6. 



3.5. Fourier Sine/Cosine Series. Finally, we consider random (truncated) 
Fourier sine/cosine series: 


fj( z ) 


cos(|z), 

sin(E^z), 


j even. 
j odd. 
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2 10 

Fig 4. Random degree 10 Taylor polynomials: rjo + ifr z + 772 fr + • • - + »?iofor- 2+ e empirical 
distribution on the right was generated using 500,000 random sums. 


The n = 2 case is shown in Figure 7 and the n = 10 case is shown in Figure 
8 . For this example, it is easy to compute the key functions. Assuming that 
n is even and letting m = n/2, we get 


m 


Aq(z) = 771 + 1 

Bq(z) = m + 1 + 2 sinh 2 (j?/) 
j= 1 

Ai {z) = 0 

m 

Bi (z) = -2ij^j cosh (jy) sinh (jy) 

j= 1 

A 2 (z) = 777.(771 + 1) (2m + l)/6 

m 

B 2 (z) = 2 ^ j 2 sinh 2 (jy) 


3= 1 


From these explicit formulas, it is easy to check that the density function 
h n (z ) depends only on the imaginary part of z as in evident in Figures 7 and 
8 . It is also easy to check that the distribution on the real axis is uniform; 
i.e., the density function g n (x) is a constant: 

9n(x) = ^yJn(n + l)/3. 

4. Final Comments and Suggested Future Research.. The ma¬ 
chinery developed in this paper can be applied in many situations that we 
have not covered. For example, if the coefficients are assumed to be inde¬ 
pendent complex Gaussains (instead of real), then we can apply the same 
methods and we expect that the computations will be simpler. In this case, 
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- 2.5 2.5 


Fig 5. Random degree 10 root-binomial polynomials: rj o + Vi\f^x z + *72 \J + 
r/s yj z 3 + ■ ■ ■ + r/ioz 10 . The empirical distribution on the right was generated using 

50,000 random sums. 


the intensity function does not have mass concentrated on the real axis (i.e., 
g n = 0) and the intensity function is rotationally invariant. 

Acknowledgement. The author would like to thank John P. D’Angelo 
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APPENDIX A: ALGEBRAIC SIMPLIFICATION OF THE FORMULA 

FOR dF/dZ. 

Substituting the derivatives given in (21) and (22) into the formula (20) 

for dF/dz, we get that the denominator simplifies to 

= (.Bo-Do + Bq — AqAoJ 
= (^BqDq + D^j 


denom ( 


(B 0 + Do) 2 D 2 0 . 
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Fig 8 . Random sum of first 11 terms in a sine/cosine Fourier series: r/o + 771 sin( 2 ) + 
r /2 cos(z) + • • ■ + 779 sin(5z) + 1710 cos(52). The empirical distribution on the right was gen¬ 
erated using 26,000 random sums. 

and the numerator of the formula becomes 

num {jppj = ( B o D o + B o - AqAq) (^B 2 D 0 + B 1 B o B i ^ A o A i 

+B\B\ + BqB 2 — 2A\Ai^j 

—{B\Dq + B 0 B 1 - AoAipBiDo + Bo B o B i~ A oAi 

\ D 0 

-\-2BqB\ — 2A\Aq^ . 

The first step to simplifying the numerator is to replace Bq — AqAq in the 
first term with Dq (like we did in the denominator) and factor out a 1/Dq 
to get 

num = D o)D 0 (b 2 D 2 0 + B\BqB\ — B\AqA\ 

+B\BiDq + BqB 2 Dq — 2AiA\Dq S j 
— — ((Bo + Dq)Bi — AqAA^D 2 + BqBqB\ — B\AqA\ 
+2BqB\Dq — 2A\AqDq S j . 
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Next, we bundle together the terms that have a5o + Do factor: 

( dF \ 1 

= + Dq)\B2Dq + BiBqB\Dq — BiAqA\Dq 

+B\BiDq + BqB 2 Dq — 2A\AiDq 

~\Bi\ 2 Dq - DoDolDil 2 + B 2 A 0 A 1 

— 2Do|Di| 2 Do + 2A\AqBiD^j 

+—^4o^4i (B\Dq + BqBqBi — BiAqA\ 
d o v 

+ 2 D 0 -B 1 D 0 — 2A\AqDq s j. 

Now, there are several places where we can find Do+Do factors. For example, 
the big factor containing eleven terms can be rewritten as follows: 

B 2 D$ + D 1 D 0 D 1 D 0 — Di^lo^iDo + B\B\Dq + BqB 2 Dq — 2A\A\Dq 
— |Di| 2 Dq — DoDo|Di| 2 + B 2 { AqA\ — 2Do|Di| 2 Do + 2^4iv4oDiDo 
= (D 0 + D 0 )(D 2 D 2 - D 0 |Di| 2 + AqA\B\) - 2|^| 2 D 2 . 

We also look for Do + Do factors in the five-term factor: 

DiDq + D 0 D 0 D 1 — B\AqA\ + 2 D 0 D 1 D 0 — 2A\AqDq 

= (Do + Do) ((Do + Do)Di - AqA\) — AqA\Dq. 

Substituting these expressions into the formula for the numerator, we get 

D 0 num = (D 0 + D 0 ) ((D 0 + D 0 )(D 2 D 2 - DqID^ 2 + AqA^) - 2\A^D 2 

+A 0 Ai ((D 0 + D 0 ) ((D 0 + D 0 )Di - AqA\) - AqA\Dq) . 
Rearranging the terms, we see that 

( dF\ 

vjwJ = (Do + Do) 2 (D 2 Dq — Do|Di| 2 + ylo^iDi + Ro^iDi) 

-(Do + Do) (2|A 1 | 2 D 0 2 + |A 0 | 2 |^ 1 | 2 ) 

— |^4o| 2 |^-i| 2 D 0 . 



22 


ROBERT J. VANDERBEI 


Here’s the tricky part... replace Dp in the second row with Bp — |Ho| 2 and 
the second and third row simplify nicely: 

(B 0 + D 0 ) (2|H 1 | 2 j D 2 + IHolVxl 2 ) + |^4 0 | 2 |^ 1 | 2 £>o 

= (Bo + Do) (2|H 1 | 2 H 2 - |A)| 2 |Hi| 2 ) + \Ao\ 2 \A^D 0 
= 2|H 1 | 2 D 0 3 - |H 0 | 2 |H 1 | 2 Ho + 2\A\\ 2 BqDq 
= \AtfBo (2 BI-\Ao\ 2 + 2BoDo) 

= \A x \ 2 Bo (Bo + Dof. 

Substituting this expression into our formula for the numerator, we now 
have 

L> 0 num = (Bo + Do) 2 (B 2 D 2 -B 0 (\A 1 \ 2 + \B 1 \ 2 ) + AoA 1 B 1 + AoA 1 B 1 ). 

Finally, we get a simple formula for dF/dz: 

8F = B 2 D 2 - BodA ^ 2 + |Di| 2 ) + ApArBr + ApA l B l 
dz Djj 
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